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ABSTRACT 

Three pairs of possibly correlated ultra-high energy cosmic ray events were 
reported by Hayashida et al (1996). Three different numerical models are 
combined to study the propagation of the corresponding particles through both 
the intergalactic and galactic magnetic fields. The spatial dependences of fields 
and galaxies are accounted for. The results suggests that the pairs are chance 
clusterings. 

Subject headings: Cosmic Rays — large-scale structure — magnetic 
fields 
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1. Introduction 

Recently, Hayashida et al. (1996) reported the possible clustering of some of the ultra 
high energy events of the AG AS A experiment. If these ultra high energy cosmic rays 
(UHECR) are charged particles, protons as it is more likely, then these pairs impose severe 
constraints on the characteristics of the propagation region and/or their sources (e.g., 
Cronin 1996, Sigl et al 1996). Catastrophic extragalactic events, hke 7-ray bursts (GRB) or 
the decay of topological defects (TD), which arc able to produce the particles over a very 
short period of time, should only be consistent with the data for a suitable combination 
of low intergalactic magnetic field (IGMF) and distance to the source. Nevertheless, the 
stirring of the intergalactic medium (IGM) by large agglomerates of galaxies, shocks excited 
in binary collisions of galaxies or the bow shocks preceding fast moving galaxies in dense 
IGM environments are examples of quiescent sources that could produce chance pairings 
of UHECR events on the sky. If these quiescent sources are traced by the distribution 
of luminous matter in the nearby universe, which is known, then the probability of the 
corresponding chance pairing can be estimated and compared with the observations. 

In this letter, the results of three different calculations are presented. First, the 
trajectories of the individual particles through the galactic magnetic field (GMF) are 
calculated for each pair under different assumptions for the GMF (Medina Tanco et al, 
1997a, Medina Tanco 1997a). In the case of catastrophic events (i.e, almost simultaneous 
particle emission) this constrains the amount of time delay due to intergalactic propagation 
alone and, consequently, the range of IGMF values and source distances allowed. The 
separation angle between the momenta of the particles at their arrival at the border 
of the halo, Ohaloj can also be estimated. This is a matching condition that must be 
satisfied by the particle trajectories at the border of the halo. Second, the same numerical 
scheme of Medina Tanco et al (1997b) is used to estimate the arrival relative-deflection 
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distribution function for some allowed combinations of IGMF and distance to the source. 
The comparison of this distribution function with the previously calculated Ohalo, gives 
a quantitative idea of the likelihood of the observed events being the result of point-like 
catastrophic sources. Third, the actual distribution of extragalactic objects, as given by 
the CfA catalogue (Huchra et al 1995), is assumed to track the UHECR sources and 
to modulate the intensity of the IGMF. Consequently, with the aid of numerical three 
dimensional simulations, an all-sky arrival probability distribution function of UHECR is 
built (Medina Tanco 1997b,c) and compared with the observations. 

2. Numerical models and discussion of results 

Three different codes are used in the present work. The first one allows the calculation 
of the trajectory of an UHECR particle of known mass and charge between the border of the 
galactic halo and the detector at Earth. A complete description of the model can be found 
in Medina Tanco et al (1997a) and Medina Tanco (1997a). The results depend, of course, 
on the model used to describe the large scale galactic magnetic field. This is certainly a 
largely unexplored area. However, we expect that a rough description, satisfactory for the 
present treatment, can be attained by the models of Stanev (1997). We adopt the same 
two extreme combinations of Stanev (1997) (see also Sofue et al. 1986, Beck et al, 1996): 
(1) a bisymmctric GMF model with field reversals and odd parity (BSS-A) and (2) an 
axisymmetric GMF model without reversals and with even parity (ASS-S). The effects of a 
small Bz — 0.1 //G component are also studied in each case. 

Table 1, adapted from Hayashida et al. (1996), lists the proposed clusters of events. 
Atarr, IS the arrival time delay. The pairs were classified as type A and B, according to 
the arrival order of the highest energy particle. Only type A events, where the highest 
energy particle arrives first, can be originated in a bursting source in which the particles are 
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simultaneously released. In type B events either the source is quiescent or there is a finite 
acceleration time involved that delays the production of the high energy particle. 

Table 2 summarizes the results for pair 1 under different GMF configurations. For 
an almost simultaneous release of the particles at the sources, all the GMF configurations 
but one, imply that the source of pair 1 should lie inside the galactic halo. The maximum 
distance to the source for each one of these GMF models is indicated in column 3. Only 
the ASS-S model without a component allows an extragalactic (EG) source. In the later 
case, a maximum arrival time delay AtioM ^ 0.6 yr is left for the intergalactic portion of 
the trajectories of both particles. 

The arrival time delay between a proton of energy E and a photon can be estimated as: 

~ 9 X 10^ X f — X ( — X ( -^y' X ( H (1) 
^ V10-9G; \3QMpcJ KlO^OeVj \lMpcJ ^ ^ ^ 

(c.f. , Waxman and Coppi, 1996), where B, Lc and D are the intensity of the IGMF, its 
correlation length and the distance to the source respectively. If the correlation length is 
known, Eq. (1) can be used to estimate maximum IGMF for a given D and Atarr between 
two protons. Two fiducial distances have been selected for quantification purposes: D = 3 
and D = 30Mpc. The maximum values of the IGMF for these distances are listed in 
Table 2 for — 1 Mpc (Kronberg 1994, 1996). These are the constraints set upon the 
intergalactic propagation region and UHECR bursting sources by the observed pair 1, after 
considering the propagation of the particles through GMF. However, another constraint 
must be satisfied: the angle between the momenta of the particles arriving at the border of 
the halo from the IGM should be equal to the calculated 9halo in Table 2. 

To this end, numerical simulations (Medina Tanco et al 1997b) were carried out 
emulating the intergalactic propagation of the components of pair 1. Lc ~ 1 Mpc is 
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assumed, while IGMF values and distances to the sources are those of Table 2. Protons 
are injected at the sources with an E~'^ spectrum and the energy losses included are 
redshift, pair production and photo- pion production (Berezinsky and Grigor'eva 1988). The 
resulting distribution functions for the relative time delay and the arrival angle between 
the proton components of pair 1 are shown in figures 1 and 2. The average time delay 
between both protons, as given by the simulations (figure 1), is consistent with equation (1), 
although there is a considerable dispersion. Furthermore, figure 2 shows that a Oralo — 2° 
separation, as inferred for pair 1 at the border of the halo, is at the wing of the distribution. 
Therefore, if a bursting source were responsible for this pair, a very low probability event 
was observed indeed. 

Pairs 2 and 3 are type B events. This means that a point source cannot have emitted 
the UHECR simultaneously. Therefore, if the point source hypothesis is to be maintained, 
we must assume either that the source is quiescent or, if bursting, that a finite acceleration 
time is involved which delays the emission of the high energy component. In this case, the 
sum of the arrival time delay and the time delay due to propagation through the GMF and 
IGMF, is a lower limit to r^, the lifetime of the source. Again, the galactic and intergalactic 
trajectories must verify the matching of 9halo a-t the galactic halo border. It is found that 
^_ffyiLo(pair 2) ~ 2° and ^^_f/ALo(pair 3) ~ 2° — 5.5", depending on the GMF model adopted. 
Numerical simulations for the IGM propagation of the proton components of pairs 2 and 3 
are also shown in figures 1 and 2. D = 30 Mpc, and Bjgm — 10~^^ and 10~^ Gauss, were 
used. The lower value of the IGMF is the one imposed by a bursting pair 1, and the second 
is the current upper limit for the IGMF. It can be seen from figures 1 and 2 that, as for 
pair 1, a 10~^^ Gauss IGMF leads to a very low probability for an event with Oralo on the 
order of a few degrees. Taking into account the galactic propagation, the lower limits for 
the hfetime of single sources for pairs 2 and 3, with Bigm ^ 10~^^ Gauss, are ~ 10 and 
~ 100 yr respectively. On the other hand, from the point of view of 9ralo, a consistent 
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picture can be obtained for a higher value of the IGMF, say near 10~^ Gauss. However, 
Ts > 10^ yr and so a single source should be quiescent and, probably, extended perhaps 
enclosing more than one galaxy in order to confine ~ 10^'^ eV particles. 

The previous results seem to point to a chance clustering of the three pairs of events, 
despite the chance probability for the pairs quoted by Hayashida et al. (1996) being only 
2.9%. We note, however, that this chance probabihty was derived under the assumption that 
the arrival direction distribution is uniform over the sky. This is arguable. Several classes of 
potential extragalactic sources have been proposed (e.g., Kewley et al 1996, Prothcroc and 
Johnson 1996, Biermann et al 1996, Halzen 1997), and these are not uniformly distributed 
over the sky. The inhomogeneity of the source's distribution should be more noticeable 
because the interaction of UHECR with the cosmic microwave background (CMB) imposes 
an upper limit D^nax ~ 10^ Mpc. Even if the actual sources are unknown, we can naively 
assume that they follow the distribution of galaxies (i.e., luminous matter) in the nearby 
universe. This is compatible with isolated galaxies, interacting galaxies, galactic bowshocks 
in high density IGMs and extended sources in turbulent IGMs powered by concentrations 
of galaxies. 

Except for some few observational determinations and upper limits (e.g., Arp 1988, Kim 
et al. 1989, Kronberg 1994) or numerical simulations of cosmological structure formation 
(Biermann 1996 and references there in) we know very little about the IGMF. These 
constraints, however, point to an IGMF structure that follows the distribution of matter 
(galaxies). Therefore, a high degree of inhomogeneity can be expected, with relatively high 
values of Bjqm over small regions (~ 1 Mpc) of high matter density (c.f., Arp 1988, Kim et 
al. 1989), pervading vast low density/low Bjom regions with Bjqm < 10~^ G. 

Following Medina Tanco et al (1997b), it is assumed that the UHECR are protons, 
and that their sources are extragalactic and hosted by, or associated with, normal galaxies. 
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It is further assumed that the magnetic field scales as n^^; , where Ugai is the local density 
of galaxies as derived from the CfA redshift catalogue (Huchra et al 1995). The IGMF is 
considered as organised in cells of size Lc of homogeneous field, such that the orientation of 
BiGM between adjacent cells is uncorrelated. Lc relates to the IGMF through the expression: 
Lc(r) cx [BiGM{f')]~'^, and the normahsation condition Lc ~ 1 Mpc for Bigm ~ 10~^ G is 
adopted. UHECR protons are injected at the galaxies with an energy spectrum oc E~'^ and 
propagated non-diffusively through the above scenario, while loosing energy via redshift, 
pair production and photomeson production (Berezinsky and Grigor'eva 1988). 

The results are displayed in the form of all-sky UHECR images of the celestial sphere 
for galaxies located at 20 < D < 50 Mpc (Figure 3.a) and 50 < i:> < 200 Mpc (Figure 3.b) 
for arriving protons with E > A x 10^^ eV. These surfaces should be representative of the 
arrival probability of UHECR at the Earth position in the Galaxy. The curved lines bound 
the region of the sky where AGASA is beheved to be sensitive (Uchihori 1996). 

We can see that the arrival probability is by no means isotropic. Furthermore, pair 2 is 
on top of a maximum of the arrival probability for sources located between 20 and 50 Mpc, 
while pair 1 is also located on a high arrival probability region for sources at more than 50 
Mpc. This is in contrast with the chance probability estimated by Hayashida et al. (1996), 
and points to either different uncorrelated sources of the components of each pair, or to very 
extended quiescent sources involving several galaxies. We also note that the sensitivity of 
AGASA is rather low in the vicinity of pairs 1 and 2. Therefore, an instrument with more 
uniform coverage (like the proposed Auger project) should probably detect an extended 
region of excess UHECR fiux at the position of the pairs. 

The third pair comes from a region of space where no large clustering of galaxies exist 
up to the depths considered. As the components cannot have originated simultaneously at 
the same extragalactic source because of galactic propagation constraints, they must have 
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come from isolated sources. This seems to indicate that very large agglomerates of galaxies, 
large enough to give a signature in figure 3, are not needed in order to accelerate UHECR. 

3. Conclusions 

The constraints deduced from the propagation of the components of the pairs of 
UHECR events proposed by Hayashida et al (1996) through the galactic and intergalactic 
medium have been analysed. 

In the case of pair 1, the low value of the IGMF, imposed by the arrival time delay 
between the protons, is inconsistent with the deflection angle between the momenta of the 
particles at the border of the halo, inferred from their galactic propagation. This makes a 
single, bursting source very unlikely. 

If the components of pairs 2 and 3 originate in common sources, then the lifetimes of 
the sources are probably larger than few times 10^ yr and, therefore, extended. This picture 
is consistent with an IGMF value not much smaller than the presently accepted upper 
limit of 10~^ Gauss (Kronberg 1996) and a distance to the sources of ~ 30 Mpc. In fact, 
the actual distribution of galaxies (Huchra et al, 1995) presents a local maximum at about 
that distance in the direction of pair 2. Furthermore, our simulations point to a maximum 
in the arrival distribution of UHECR at exactly the same position of pair 2, when sources 
between 20 and 50 Mpc are considered. Pair 1 is also located inside a maximum of the 
arrival distribution, favouring chance pairing between the components. 

This work was done with the partial support from the Brazilian agency FAPESP. 
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Figure Captions 

Figure 1: Distribution function of arrival time delays between the observed pair of 
protons in clusters 1, 2 and 3 due to propagation in the IGMF alone (i.e., at the external 
border of the galactic), halo). The simulations for pair 1 correspond to the two fiducial 
scenarios of Table 2, and match the constrain in time delay imposed by the galactic portion 
of the tracks: Bjgm = 10"^^ G and D = 3 Mpc (dotted line) and Bigm = lO'^^ G and 
D — SO Mpc (continuous hue). For pairs 2 (broken-dotted hue) and 3 (broken hue) two 
possible scenarios are explored: D — SO Mpc and Bjqm — 10~^^ G and Bjqm — 10~^ G. 

Figure 2: Distribution function of the angle between the momenta of the observed 
particles in each pair, at their arrival at the external border of the galactic halo after 
propagation through the IGMF. The conditions are the same as in figure 1. Also indicated 
is a separation angle of 2° typically obtained from the calculations of galactic propagation 
for all the three pairs. 

Figure 3: Arrival distribution of UHECR simulated under the assumption that the 
luminous matter in the nearby Universe tracks the distribution of the sources of UHECR, 
and modulates the intensity of the IGMF (see text for details). Redshift, pair production 
and photo-pion production losses are included. The arrival distributions due to sources at 
two different depths are shown: (a) < D < 50 Mpc, (b) 50 < D < 200 Mpc. Pairs 1 and 
2 lie on top of regions of high arrival probability strengthening the possibility of chance 
pairing. The solid lines bound the region of the sky actually seen by AGASA. Sensitivity is 
poor near these lines 
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Table 1: Possible Clusters of UHECRs observed by AG ASA experiment (adapted from 
Hayashida ct al. 1996) 



Pair No. 


Date 




Energy [eV] 


Type 


Igal 


bgal 


1 


93/12/03 


1.90 


210 


A 


131.2 


-41.1 




95/10/29 




51 




130.2 


-42.3 


2 


92/08/01 


2.49 


55 


B 


143.5 


56.9 




95/01/26 




78 




145.8 


55.3 


3 


91/04/20 


3.21 


43 


B 


77.9 


18.6 




94/07/06 




110 




77.6 


21.1 



Table 2: Pairl: Constraints from galactic and intergalactic propagation. 



Bgal 




Dmax [kpc] 


At/GM [yr] etHALO [1 




ASS-S 


= 


EG 


0.6 2 


IQ-ll lQ-12 






7 


4.5 




BSS-A 


= 


8 


2 






7^0 


13 


0.5 
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